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Abstract 

In this work, the analytical solution of the hyper-radial Schrodinger equation for the spherical 
Woods-Saxon potential in D dimensions is presented. In our calculations, we have applied the 
Nikiforov-Uvarov method by using the Pekeris approximation to the centrifugal potential for arbi- 
trary / states. The bound state energy eigenvalues and corresponding eigenfunctions are obtained 
for various values of n and / quantum numbers. 
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I. INTRODUCTION 



An analytical solution of the radial Schrodinger equation is of high importance in nonrel- 
ativistic quantum mechanics, because the wave function contains all necessary information 
for full description of a quantum system. There are only few potentials for which the ra- 
dial Schrodinger equation can be solved explicitly for all n and I quantum numbers. So 
far, many methods were developed, such as supersymmetry (SUSY) [1,2] and the Pekeris 
approximation [3-8], to solve the radial Schrodinger equation exactly or quasi-exactly for 
I 7^ within these potentials. Levai and Williams suggested a simple method for construct- 
ing potentials for which the Schrodinger equation can be solved exactly in terms of special 
functions [9] and showed relationship between the introduced formalism and supersymmetric 
quantum mechanics [1]. Amore et al. presented a new method [10] for the solution of the 
Schrodinger equation applicable to problems of nonperturbative nature. In addition, they 
applied the method to the quantum anharmonic oscillator and found energy eigenvalues and 
wave functions, even for strong couplings. 

The three-dimensional radial Schrodinger equation for the sperical Woods-Saxon potential 
[11] can not be solved analytically for / 7^ states, because of the centrifugal term ~ r~^. 
However, Fliigge gave an exact expression for the wave function, but a graphical method 
was suggested for the energy eigenvalues at / = [4]. It is well known that the Woods- 
Saxon potential is one of the important short-range potentials in physics. Furthermore, 
this potential were applied to numerous problems, in nuclear and particle physics, atomic 
physics, condensed matter, and chemical physics. 

Recently, an alternative method known as the Nikiforov-Uvarov (NU) method [12] was 
proposed for solving the Schrodinger equation. Therefore, it would be interesting and im- 
portant to solve the nonrelativistic radial Schrodinger equation for Woods- Saxon potential 
for / 7^ 0, since it has been extensively used to describe the bound and continuum states of 
the interacting systems. Thus, one needs to obtain the energy eigenvalues and correspond- 
ing eigenfunctions of the one particle problem within this potential. The NU method was 
used by C. Berkdemir et al. [13] to solve the radial Schrodinger equation for the general- 
ized Woods-Saxon potential for / = 0. However, it this work, the authors made errors in 
application of the NU method, which led to no-correct results [14]. In Refs.[7,8] the ra- 
dial Schrodinger equation with the spherical Woods-Saxon potential have been solved using 
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the Nikiforov-Uvarov method and Pekeris approximation for arbitrary / states. The same 
method has been employed to solve the Klein-Gordon equation with the spherical Woods- 
Saxon potential for / 7^ states Ref.[15]. 

In this work, we solve the D dimensional hyper-radial Schrodingcr equations with spher- 
ical Woods-Saxon potential for arbitrary / states. Here the NU method [12] and the Pekeris 
approximattion is applied to find energy eigenvalues and corresponding eigenfunctions of 
the considering problem. 



II. NIKIFOROV-UVAROV METHOD 



The Nikiforov-Uvarov (NU) method is based on the solutions of general second-order 
linear equations with special orthogonal functions. It has been extensively used to solve 
the nonrelativistic Schrodinger equation and other Schrodinger-like equations. The one- 
dimensional Schrodinger equation or similar second-order differential equations can be writ- 
ten using NU method in the following form: 

+ ^^'^^^ + ^^^^^ = °' ^^-'^ 
where 17(2;) and ^{z) are polynomials, at most second-degree, and t{z) is a first-degree 
polynomial. 

Using Eq.(2.1) the transformation 

iIj{z) = ^z)y{z) (2.2) 
one reduces it to the hypergeometric-type equation 

a{z)y" + T{z)y' + Xy = 0. (2.3) 



The function ^{z) is defined as the logarithmic derivative [12] 

^'(z) 7t(z) 



(2.4) 



$(z) a{z)' 
where 7r{z) is at most the first-degree polynomial. 

The another part of ip{z), namely y{z), is the hypergeometric-type function, that for fixed 
n is given by the Rodriguez relation: 



where 5„ is the normahzation constant and the weight function p{z) must satisfy the con- 
dition [12] 

^^{a{z)p{z))=T{z)p{z), (2.6) 

with t{z) — t{z) + 27r{z). 

For accomphshment of the conditions imposed on function p{z), the classical orthogonal 
polynomials, it is necessary, that polynomial t{z) becomes equal to zero in some point of an 
interval (a, b) and derivative of this polynomial for this interval at a{z) > will be negative, 
i.e. r'{z) < 0. 

The function 7r{z) and the parameter A required for this method are defined as follows: 



/ ~ / / / ~\ 2 
a — T , I a' — T 



7r(z) = ^- ± y (^^- j -a + ka, (2.7) 

X = k + T:'{z). (2.8) 

On the other hand, in order to find the value of the expression under the square root 
must be the square of a polynomial. This is possible only if its discriminant is zero. Thus, 
the new eigenvalue equation for the Schrodinger equation becomes [12] 

A = A„ = -nr' - I^il^^", (^ = 0, 1, 2, ...). (2.9) 

After the comparison of Eq.(2.8) with Eq.(2.9), we obtain the energy eigenvalues. 

III. SOLUTIONS OF THE SCHRODINGER EQUATION WITH THE WOODS- 
SAXON POTENTIAL 

Woods and Saxon introduced a potential to study elastic scattering of 20 MeV protons 
by a heavy nuclei [11]. The spherical Woods-Saxon potential that was used as a major part 
of nuclear shell model, has received a lot of attention in nuclear mean field model. The 
spherical standard Woods-Saxon potential [11] is defined by 

no = -7—^7^^' « « (3.1) 
1 exp y—^) 

This potential was used for description of interaction of a neutron with a heavy nucleus. The 
parameter Rq is interpreted as radius of a nucleus, the parameter a characterizes thickness 
of the superficial layer inside, which the potential falls from value V — Q outside of a nucleus 
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up to value V — —Vq inside a nucleus. At a = 0, one gets the simple potential well with 
jump of potential on the surface of a nucleus. 

Using D-dimensional {D > 2) polar coordinates with polar variable r (hyperradius) anf 
angular variables 61,62, ■ ■ ■ ,^z)-25 (hyperangles), the Laplasian operator in polar coordi- 
nates (r, 6*1, 6*2, ... , 9d-2-i 4>) of is 



^ dr \ dr J ' 

where A|, is a partial differential operator on the unit sphere S^~^ (Laplace-Beltrami oper- 
ator, or grand orbital operator, or hyperangular momentum operator) defined analogously 
to a three-dimensional angular momentum, Avery [16]. 

The D dimensional Schrodinger equation with spherically symmetric potential V{r) has 
the form [16] 

+ V{r) - ipnim{r,nD) = 0, (3.2) 
where is the reduced mass, h is the Planck's constant and 

i^nim{r, ^d) = Rni{r)Yi^{nD)- (3.3) 

The Laplasian operator divides into a hyper-radial part (^^~^^) angular 

part — i-G. 

Vl = r'-^l (r''-'l) - (3.4) 



dr \ dr J hPr"^ ' 

where L^, is the grand orbital angular momentum operator. The eigenfunctions of are 
the hyper-spherical harmonics Yimifio) 

LlYi^i^D) = hH{l + D- 2)Yi^{nn), (3.5) 

where I is the angular momentum quantum number. 

After substituting the Eqs.(3.3), (3.4), (3.5) into (3.2) and using the fact that ipnim{r, ^d) 
is the eigenfunction of L|, with eigenvalue h^lil + D — 2), we obtain an equation known as 
the hyper-radial Schrodinger equation with Woods-Saxon potential 



(PRnijr) ^ D-ldRni{r) ^ 2fi 
dr"^ r dr 

Introducing a new function 



h%l + D-2) 

El -\- 



1 + exp (^) 2//r2 



i?„KO = 0,(0<r<oo). 

(3.6) 



D-l 



Uni{r) = r 2 Rni{r), 



Eq.(3.6) reduces to 



(funiir-) 2/i 



E + 



1 + exp (^) 



Uni{r) = 0. 



(3.7) 



introducing a new parametr 

Eq.(3.7) takes the form 

d?Uni{r) 2(1 



^1 + 



D-3 



E + 



Vo 



hH (l + 1) 



1 + exp (^) 2//r2 



Uni{r) = 0. 



(3.8) 



Equation (3.8) has the same form as the equation for a particle in one dimension, except for 
two important differences. First, there is a repulsive effective potential proportional to the 
eigenvalue of fi^l{l + 1). Second, the radial function must satisfy the boundary conditions 
m(0) = and u{oo) = 0. 

It is sometimes convenient to define in Eq.(3.8) the effective potential in the form: 

hH{T+ 1) 



Veffir) = V{r) + 



(3.9) 



Then, the radial Schrodinger equation given by Eq.(3.8) takes the form 

d'^Uniir) 2/1 



dr 



2 . ^[E-Veff{r)]Unl{r)^0. 



(3.10) 



If in Eq.(3.1) introduce the notations 

_ r - -Rq _ Ro 

X — , Q! — , 

Kq a 

then the Woods-Saxon potential is given by the expression 



Vws — 



1 + exp(Q;a;) 

The effective potential together with the WS potential for i 7^ can be written as 

K„W^^W + W.)^5|±ii-^^^. (3.11) 

It is known that the Schrodinger equation cannot be solved exactly for this potential at 
the value / 7^ using the standard methods as SUSY and NU. Prom Eq.(3.11) it is seen 
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that the effective potential is a combination of the exponential and inverse square potentials, 
which cannot be solved analytically. Therefore, in order to solve this problem we can take 
the most widely used and convenient for our purposes Pekeris approximation [3-8, ]. This 
approximation is based on the expansion of the centrifugal barrier in a series of exponentials 
depending on the internuclear distance, taking into account terms up to second order, so 
that the effective /-dependent potential preserves the original form. It should be pointed 
out, however, that this approximation is valid only for low vibrational energy states. By 
changing the coordinates x — ^^^^ or r — Ro{l + x), the centrifugal potential is expanded 
in the Taylor series around the point x = {r — Rq) 

V,(r) = = "^'(' + 1) , 1 = ^ (1 - 2. + 3.^ - 4.' + . . .) , (3.12) 

where ?= 

According to the Pekeris approximation, we shall replace potential Vj{r) with expression 

V{{r) = 6(co + + ^"j . (3.13) 

' V 1 + expax (l + expcKx) V 

In order to define the parameters Co, Ci and C2, we also expand this potential in the 
Taylor series around the point x — {r — Rq): 



Co + ^ + - -{Ci + C2)x + —C^x^ + 7^ (C^i + C2) x^ - -^C2x'^ + 



2 4 y 4 ' ' 16 48 ' ' 96 

(3.14) 



Comparing equal powers of x Eqs.(3.12) and (3.14), we obtain the constants Co, Ci and C2: 

4 12 8 48 48 

<-^0 — J- 1 ^! ^1 — 0' 2 — 

a a Q; Q;- a- 

Now, the effective potential after Pekeris approximation becomes equal to 

l + exp(Q;a;) (1 exp(Q;a:)) 
Instead of solving the radial Schrodinger equation for the effective Woods-Saxon potential 
Veffir) given by Eq.(3.11), we now solve the radial Schrodinger equation for the new effective 
potential given by Eq.(3.15) obtained using the Pekeris approximation. Having 

inserted this new effective potential into Eq.(3.10), we obtain 



l + exp(^) (i + exp(rz^))^ 
7 



Ur,i{r) = 0. (3.16) 



We use the following dimensionless notations: 

with real e > for bound states; /5 and 7 are real and positive. 
If we rewrite Eq.(3.16) by using a new variable of the form 



1 + exp 



we obtain 



u"{z) + -7-—y{z) + .J u{z) = 0,{0<z< 1), (3.18) 

with t{z) = 1 - 2z; a{z) = z{l - z);a{z) = -e^ + /3'^z - i^z^. 
In the NU-method, the new function 7r(z) is 

The constant parameter k can be found employing the condition that the expression 
under the square root has a double zero, i.e., its discriminant is equal to zero. Hence, there 
are two possible functions for each k: 

( (e- Je^ - /32 + 72) z-€, for k = 13^- 2e^ + 2€Je''-(3^ + -f^, 

Tr(z) = ± < > , / , (3.20) 

[U+ Ve' + 7') z-e, for k ^ /3'' - 2e^ - 2ey/e^ - + -f^. 

According to the NU-method, from the four possible forms of the polynomial n{z) we select 
the one for which the function t{z) has the negative derivative and root lies in the interval 
(0,1). Therefore, the appropriate functions n{z) and t{z) have the following forms: 

Tr{z) = e - + - + 7^) (3.21) 

t{z) = 1 + 2e - 2 (^1 + e + ^/e^ - /^^ + 72^ z, (3.22) 

and 

k^/3^-2e^- 2ey/e^ - + 72. (3.23) 
Then, the constant X = k + 7r'{z) is written as 

A = - 2e2 - 2e - ^2 + ^2 _ ^ + ^^2 _ ^2 + ^2_ (3 34) 



An alternative definition of A„ (Eq.(2.9)) is 

A = A„ = 2 (^e + - + 7') « + n{n + 1). (3.25) 
Having compared Eqs.(3.19) and (3.20) 

p^-2e^-2ey/e^ - + j^-e-^/e^ - /ja + ^2 = 2(^6+ ^ -(3^ + 72) n+n{n+l), (3.26) 
we obtain 

2 2 



e + Ve' - + 7' + ^ + ^ - J + = 0, (3.27) 



or 



e + - /32 + 72 = n'. (3.28) 
After some elementary calculus, one sees that 

/ 0^ - 7^ 

e- Ve2-/32 + y = ^^^. (3.29) 

Here 

„- = _„+VlW-_l^ (3.30) 

2 

n being the radial quantum number (n = 0, 1, 2, . . .). From Eq.(3.28) and (3.29), we find 

1 / , , /?'-7^ 



and 

Ve' -P' + l' = l{n'- (3.32) 
Because for the bound states e > and Eqs.(3.28, 3.29, 3.31, 3.32) we get 

n' > (3.33) 

and 

< - 7^ < n'l (3.34) 

If n' > 0, there exist bound states, otherwise, there are no bound states at all. By using 
Eq.(3.30) this relation can be recast into the form 



< n < — , (3.35) 



i.e. it gives the finite coupling value. 



If — 7^ > 0, there exists bound states; otherwise there are no bound states. Inequahty, 
which obtained after substituting /3, 7, Ci, C2 into Eq.(3.34), gives the definite couphng value 
for the potential depth Vq: 



hH{l + l)a 



(3.36) 



After substituting 7, 5, C2 into Eq.(3.36), we find 



< n < 



192a4;(;+l) 



1 + 



- 1 



(3.37) 



The exact energy eigenvalues of the Schrodinger equation with the Woods-Saxon potential 
are derived as 



n" + - 7 



2\ 2 



(3.38) 



Substituting the values of 5, Cq, C2, n\ (3 and 7 into (3.38), one can find E^i 

^(D) _ hH{T+ 1) /. 12a^' 



e: 



nl 



2iJLa? 



16 



+ - 



Rl 



fia'^Vo _ 4:l{l+l)a'^ 



1 2 



-"0 



+ 



(3.39) 



In three-dimensions, the case when we set the parameter value D = 3 in Eq.(3.39), then 
we obtain the energy of the Woods-Saxon potential 



^nl — 



- 2n - 1 



16 



+ 



1 



4 




2 




ft2 iig 





l + «^-2n-l 



^2 



(3.40) 



which is identical to the one obtained in Rcf. [8] . 

If both conditions (3.36) and (3.37) are satisfied simultaneously, the bound states exist. 
For very large Vq the /-dependent effective potential has the same form as the potential with 
/ = 0. Whenever D = 3, the bound states of the system do not exist in / = state, because 
(3.36) and (3.37) inequahties are not satisfied. Hence, the radial Schrodinger equation for 
the standard Woods-Saxon potential with zero angular momentum has no bound state. 
Furthermore, whenever D > 3, the bound states of the system exist in / = state. Thus, 
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the energy spectrum equation (3.39) is limited, i.e. we have only the finite number of energy 
eigenvalues. 

According to Eq.(3.39), the energy eigenvalues depend on the depth of the potential Vq, 
the width of the potential Rq, surface thickness a, and the parameter D. If constraints 
imposed on n and Vq satisfied, the bound states appear. From Eq.(3.36), it is seen that the 
potential depth decreases when the parameter Rq increases and the parameter a decreases, 
and vice versa. Therefore, one can say that the bound states exist within this potential. 

In addition, we have seen that there are some restrictions on the potential parameters 
for the bound state solutions within the framework of quantum mechanics. That is, when 
the values of the parameters Vq and n satisfy the conditions (3.36) and (3.37), we ob- 
tain the bound states. We also point out that the exact results obtained for the standard 
Woods-Saxon potential may have some interesting applications for studying different quan- 
tum mechanical and nuclear scattering problems. Consequently, the found wave functions 
are physical ones. 

Now, we are going to determine the radial eigenfunctions of this potential. Having sub- 
stituted 7r{z) and a{z) into Eq.(2.4) and then solving first-order differential equation, one 
can find the finite function ^{z) in the interval [0, 1] 



$(^) = ^^(l-z)V^^^. (3.41) 

It is easy to find the second part of the wave function from the definition of weight 
function: 



p{z) = z^' (1 - ^)V^^-/5^+7^ , (3.42) 
and substituting into Rodrigues relation (2.4), we get 



'm{z) = B^z-'^ (1 - z)-'^/^^^' ^ \z^+'^ (1 - zr-"'^/^"^' 



(3.43) 



where is the normalization constant and its value is ^ [17]. Then, yn is given by the 

Jacobi polynomials: 

yn{z)^P} ^(1-2;.), 

where 

1 fl"' r n 

pi'^'^Hi - 2^) = (1 - -2"^" (1 - ^T^i ■ 

" ^ ^ n! ^ ^ dz"" I ^ ^ J 
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The corresponding Uni{z) radial wave functions are found to be 

u^(z) = C„,z' (1 - pf'-^^^'hl - 2.). (3.44) 

where C„; is a new normahzation constant determined using J^[uni{r)]^dr = 1. 

IV. CONCLUSION 

In this paper, we have analytically calculated energy eigenvalues of the bound states 
and corresponding eigenfunctions in the new exactly solvable Woods-Saxon potential. The 
energy eigenvalue expression for Woods-Saxon potentials is given by Eq.(3.39). As it should 
be expected (see Eq.(3.39)), for any given set of parameters Vo,Rq, a and D, the energy 
levels of standard Woods-Saxon potential are positive. We can conclude that our results 
are interesting not only for pure theoretical physicist but also for experimental physicist, 
because the results are exact and more general and useful to study nuclear scattering. 
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